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I. INTRODUCTION 

The Agent Based Model (ABM) studied in this paper is a generalization of the “toy-ABM” (dubbed Mark-0) 
recently introduced in [1], following previous work by the group of Belli Gatti et al. [2]. Mark-0 considers a stylized 
economy with firms and households, but no banks, no interest rates on loans and deposits, and therefore no direct 
concept of “monetary policy”. As discussed at length in [I], the original motivation of Mark-0 was mostly to illustrate 
the importance of phase diagrams and phase transitions in the context of ABMs, in particular the sensitivity of the 
state of the (artihcial) economy on a subset of parameters. Small changes in the value of these parameters were indeed 
found to induce sharp variations in aggregate output, unemployment or inflation. In other words, endogenous crises 
can occur in such economies, as the result of insignificant or anecdotal changes in the environment. This possibility 
is quite interesting in itself, and must be contrasted with more traditional economic models, such as the popular 
DSGE framework [3, 4], where the dynamics is linear and only large exogenous shocks can cause havoc. As recently 
pointed out by O. Blanchard in a very inspiring piece [5]: We in the field did think of the eeonomy as roughly linear, 
constantly subject to different shocks, constantly fluctuating, but naturally returning to its steady state over time. 

The main lesson of the crisis is that we were much closer to “dark corners” - situations in which the economy could 
badly malfunction - than we thought. 

Because they can deal with non-linearities, heterogeneities and crises, ABMs are often promoted as possible al¬ 
ternatives to the DSGE models used in central banks as guides for monetary policy [4, 6-8]. It is therefore clear 
that introducing interest rates monitored by a central bank in Mark-0 is mandatory for policy makers to develop any 
interest in the ABM research program in general and our model in particular. The aim of the present paper is to 
parsimoniously extend Mark-0 as to capture the effects of monetary policy on the course of the economy. We first 
identify and model several channels through which interest rates can feed into the behavior of firms and households. 
We then study different policy experiments, whereby the “Central Bank” attempts to reach a target inflation and/or 
unemployment level using a Taylor rule to set the interest rate (see Eq. (3) below). We find that provided the economy 
is far from phase boundaries (or “dark corners” [5]) such policies can be successful, whereas too aggressive policies 
may in fact, unwillingly, drive the economy to an unstable state, where large swings of inflation and unemployment 
occur. 

Our Agent-Based framework is voluntarily bare bones. It posits a minimal set of plausible ingredients that are 
most probably present in the real world in one form or the other. Eor example, we assume reasonable heuristic rules 
for the hiring/firing and wage policies of firms confronted with over- or under-production, or with a rising level of 
debt. Similarly, our model encodes in a schematic manner the consumption behavior of households facing inflation 
and rising rates, that is in fact similar to the standard Euler equation for consumption in general equilibrium/DSGE 
models [3]. Our approach is therefore prone to the usual critique addressed to ABMs: the rules we implement are 
- although reasonable - to some degree arbitrary. The ABM community is of course aware of this weakness, with 
many attempts to resolve it, such as imposing consistency constraints on behavioral assumptions (see [9], Appendix 
A), or, even better, relying on micro-panel data that reveal how Hrms and households actually make decisions under 
different macro-economic or specihc conditions (see for example [10, 11] for the behavior of firms, and [12-14] for the 
behaviour of households). However, reliable empirical data are still rather scarce and do not allow yet to answer all 
the questions needed to constrain and calibrate an ABM, even simple ones like Mark-0. 

Our philosophy, explained in detail in [1], is different. We argue that the qualitative, aggregate behavior elicited by 
the Mark-0 model is in fact robust and generic, although the actual quantitative aspects may not be (as, for example, 
the precise value of the parameters of the Taylor rule beyond which instabilities occur). In other words, if one forgoes 
the idea of quantitatively predicting the macro-behaviour of the economy but is satisfied (at least temporarily) with 
a qualitative description of the possible aggregate behaviour, some progress is possible without detailed knowledge of 
the micro-rules. Our belief is backed by the idea - pervasive in many areas of science - that the aggregate properties 
of interacting entities can be classified in different phases, separated by phase boundaries across which radical changes 
of the emergent behavior take place. This idea has a long history, in particular in economics - remember, for example, 
the title of Thomas Schelling’s famous book: Micromotives and Macrobehaviour [15]; for more recent progress see 
e.g. [16-21] and for recent reviews: [22, 23]. To bolster our belief that emergent aggregate properties are robust against 
micro-changes, we have tested many variants of the model presented below and indeed found that the overall behavior 
of our artificial economy is remarkably robust - in particular the presence of instabilities and crises. Eollowing up on 
O. Blanchard’s lament [5], the existence of large swaths of the parameter space where the economy is prone to violent 
crises seems to be an unavoidable fact that we have to learn to confront with [16, 24]. 

The work presented in this paper is part of a growing literature of macroeconomic agent-based model [25, 26]. In the 
recent years several such models have been developed, allowing both to reproduce macro-stylized facts and to study 
policy design [27]. Building upon [7], Dosi et al. [28, 29] characterize the impact of fiscal and monetary policies on 
macroeconomics fluctuations. Fiscal policies are found to have a greater role in dampening business cycles, reducing 
unemployment and the likelihood of economic crises. Another example somewhat related to our work is Ref. [30], 
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which focuses on a simple ABM of financial markets coupled with a New Keynesian DSGE model, which allows one 
to reproduce endogenously stock price bubbles and business cycles and to study the introduction of financial taxes. 
Calibration and validation procedures have been explored in [31-33]. 

Although similar in spirit, the present study builds upon a slightly different perspective, following the framework and 
philosophy developed in [Ij. In particular, we model an idealized closed economy with linear production capabilities, 
no capital (labor is the only input for production), no innovation and growth, no financial sector, and only a minimal 
set of additional behavioural rules which couple the Central Bank monetary policy to firms’ and agents’ choices. In 
this sense, our model is much more parsimonious than what is found in the recent ABM literature and voluntarily 
overlooks important actors and processes at play in the real world. Our central concern is instead to characterize 
the qualitative behaviour of the emerging economy and investigate its dynamical stability with respect to the Central 
Bank’s policy. 

The outline of the paper is as follows. In Section II we give a brief description of the Mark-0 model proposed 
in [I] and detail the minimal additional rules we need to couple the monetary policy with firms’ and agents’ decisions. 
In Section III, we investigate the “natural” state of our toy-economy, in the presence of interest rates but without 
any Central Bank intervention. In Section IV, we perform several policy experiments: an expansionary monetary 
shock, and the implementation of a Taylor-rule monitoring of the interest rate by the Central Bank to achieve given 
employement and inflation targets. We establish the phase diagram of our model in the presence of this Taylor-rule 
based intervention. We briefly compare our findings to the prediction of DSCE models. We conclude in Section V. 
Appendix A gives further details about Mark-0; Appendix B investigates the continuous phase transition induced by 
the wary behaviour of indebted firms; finally. Appendix C gives a full pseudo-code of our model that should allow 
easy duplication of our results. 


II. DESCRIPTION OF THE MODEL 
A. Brief summary of the minimal “Mark-0” model 

The Mark-0 closed economy is made of firms and households. While the latter sector is represented at an aggregate 
level, firms are heterogeneous and treated individually. Each firm i = 1,... ,Np at time t produces a quantity V(t) of 
perishable goods that it attempts to sell at price pi (t) . It needs a number of Ni (t) = Yi (t) / Q of employees^ to produce 
Yi{t), and pays a wage Wi(t). The demand Di for good i depends on the global consumption budget of households 
CB{t), itself determined as a fraction of the household savings (that include the last wages), and it is a decreasing 
function of the asked price Pi{t), with a price sensitivity parameter that can be tuned - see Appendix A. 

To update their production, price and wage policy, firms use reasonable “rules of thumb” [I] that we detail in 
Appendix A and that were already (partly) justified in the original work of Belli Catti et al. [2]. Eor example, 
production is decreased and employees are made redundant whenever Yi > Di, and vice-versa.^ The adjustment 
speed can however be asymmetric, i.e. the ratio R of hiring adjustement speed to firing adjustement speed is not 
necessarily equal to one, for example because of labour laws. This turns out to be one of the most important control 
parameter that determines the fate of the overall economy. 

In the initial state of the economy, firms are heterogenous in size and prices (with a uniform distribution around 
the average size and price, see Appendix C for details), but all offer the same wage and see the same demand. This 
is arbitrary and of little importance, because the stationary state of the model is statistically independent of the 
initialization. 

When the Mark-0 economy is set in motion, it soon becomes clear that some firms have to take up loans in order 
to stay in business. One therefore immediately has to add further rules for this to take place. In the zero-interest rate 
world of Mark-0, we let firms freely accumulate a total debt up to a threshold that is a multiple 0 of total payroll, 
beyond which the firm is declared bankrupt (its debt is then repaid partly by households and partly by surviving firms, 
such that there is no net creation of money). Erom this point of view the parameter 0 determines the maximum credit 
supply available to firms. Eixing the value of 0 plays the role of a primitive monetary policy, since the total amount 
of money circulating in the economy (‘broad money’) directly depends on 0 [I]. When 0 = 0, no debt is allowed 


^ i^i is the productivity of firm i. We chose = 1 in [1] and we will stick to this choice throughout the present paper as well. 

^ As a consequence of these adaptive adjustments, the economy can reach (in some regions of the parameter space) equilibrium, corre¬ 
sponding to the market clearing condition one would obtain in a fully representative agent framework. However, fluctuations around 
equilibrium persists in the limit of large system sizes giving rise to a rich phenomenology, see [1, 2]. 
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FIG. 1: (Left) Phase diagram in the R — Q plane of the basic Mark-0 model as obtained in [1] with wage update. There 
are four distinct phases separated by critical lines. The Full Employment (FE) phase {R > Rc, O large) is characterized by 
positive average inflation, while there is deflation in the Full Unemployment (FU) phase {R < Rc). Endogenous Crises (EC) 
are characterized by alternating cycles of inflation and deflation, and occur for R > Rc, 0 intermediate. Finally, R > Rc, small 
0 correspond to a region of small inflation and Residual Unemployment (RU). The location of phase boundaries is only weakly 
affected by the choice of the other parameters of Mark-0, see [1]. [Right) Typical trajectories of the unemployment rate u[t) for 
each of the phases. In the inset, the price variations are shown, displaying either inflation (in the FE phase) or deflation (in the 
EU phase). 7 j 7 ^ (resp. 7 ™^) sets the characteristic time scale for price adjustments (resp. wages) in the model, see Appendix 
A. The surprising occurrence of endogenous oscillations in the EC phase can be fully understood analytically, see [34]. 


(zero leverage), while when 0 —> 00 , firms have not limit on the loans they need to continue business (unbounded 
leverage). 

While there are several other parameters needed to define completely Mark-0 (9 parameters in total, see Appendix 
A), the detailed investigation of [1] has established that only R and 0 are the ones determining the phase-diagram of 
the model, shown in Fig. 1, where we also plot typical trajectories of the economy in each phase. It is important to 
stress that this diagram is extremely robust against both details of the model specification and the value of the other 
parameters, which only affect the above phenomenology quantitatively, but leave the qualitative emergent behavior 
essentially unchanged. Its salient features are [1]: 

• When 0 = 00 the economy is characterized by two distinct phases separated by a first order (discontinuous) phase 
transition as a function of the parameter R. When R < Rc (i.e. fast downward production adjustments), one 
finds at long times a collapse of the economy towards a deflationary/low demand/full unemployment state (FU). 
For R ^ Rc, on the other hand, the long run state of the economy is characterized by a positive inflation/high 
demand/full employment phase (FE). 

• When 0 < 00 the above description holds but must be refined to allow for the appearance of three sub-phases 
for R > Rc'. 

1. a full employment and inflationary phase for high values of 0 (the FE phase, similar to the 0 = 00 case); 

2. a phase for intermediate values of 0 characterized by high employment and inflation on average, which is, 
however, intermittently disrupted by “endogenous crises” (EC), that temporarily bring deflation and high 
unemployment spikes; 

3. a phase with zero inflation and residual unemployment for small 0 (the RU phase), where the impossibility 
to obtain loans creates a positive stationary level of bankruptcies. 

Note that the unexpected occurrence of purely endogenous oscillations in the EC phase can in fact be understood 
fully analytically [34] and it is a robust feature that only relies on very mild assumptions about the destabilizing 
feedback mechanisms present in the economy. 
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B. Introducing interest rates in Mark-0 

We now introduce in the model a banking system made up of a Central Bank (CB) which sets the base interest 
rate poit) (and, as part of a prudential policy, the parameter 0 that controls the maximum credit supply available to 
firms), and a private banking system that will act as a transmission belt for the CB policy, by setting interest rates 
on deposits (p'^(t)) and on loans {p^{t)). These interest rates will in turn impact the economy through three channels 
that we detail below: a) direct cost of loans and gains on deposits; b) behavior of the firms; c) behavior of households. 


1. The Central Bank policy 


The CB attempts to steer the economy towards a target inflation level tt* and employment level e* (equivalent to a 
target output, since the productivity C is set to unity in the present version of the model). The instantaneous inflation 
7r(t) and employment level £{t) are defined here as: 


7r(t) 


p{t - 1) 


pit) 


Y.^p^it)y^it) 

E.y^it) 


where p{t) is the production-weighted average price, and: 

^it) = ’ u(t) = l-£(t), 

I 


( 1 ) 

( 2 ) 


where N is the total workforce, and £{t),u(t) are respectively the employment and unemployment rate. 

We will assume that the monetary policy followed by the CB for fixing the base interest rate is described by a 
standard Taylor rule of the form [3, 35]:^ 

Poit) = max(/9* -t 10 - TT*] -t ^elog [£(t)/£*], 0) (3) 


where p* is the “natural” interest rate that would prevail if the target inflation tt* and target employment s* were 
reached, and aTr,s > 0 quantify the intensity of the policy (high values of the parameters correspond to aggressive 
policies). Note that po is constrained to be non-negative. The factor 10 in front of </>.„• is only there for convenience, 
such that interesting values of and (pe are of the same order of magnitude. The notation x(t) corresponds to the 
exponential moving average of the variable x{t), defined as:^ 


x{t -I- 1) = ujx{t) -I- (1 — uj)x{t). 


( 4 ) 


In order to avoid unnecessary excessive policy response when the target employment rate is too far from the actual 
employment rate, we actually define a one-time-step target employment rate e* as 

£*= min{1.025 £(t),£*} (5) 

meaning that if the employment rate is much lower than the policy target e* the CB will try to increase it by 2.5% 
at each time step until the target is reached. Such a regularization was found not to be needed for inflation. 


2. The banking seetor 

Households in the Mark-0 economy cannot borrow and are thus characterized by their total savings S{t) > 0. 
Firms, on the other hand, can have either deposits {Si > 0) or liabilities {Si < 0). Defining = ^^max(fi,0) and 
S~ = — ^ ■ min {Si, 0), the balance sheet of the banking system reads: 

M^S-{t)=S{t)^S+{t) = X{t), (6) 


® In Gall’s reference book [3], the quantity be is noted <f)y. 

^ We chose uj = 0.2, which corresponds to an averaging time of about —l/log(l — lOt) ~ 4.5 time steps. 
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where M is the amount of currency (or initial deposits) created by the central bank, which is kept fixed in time, and 
X is the total amount of deposits, therefore to initial deposits M plus the money created by the banking system when 
issuing loans. 

We now assume that the banking sector fixes the interest rates on deposits and loans and p^{t) respectively) 

uniformly for all lenders and borrowers according to the following rules:® 

P^(t) = Poit) + ■I' 
d/,^ _ Po{t)^~{t) - (1 - 

^ X{t) ■ ^ ^ 

where T>{t) is the aggregate costs coming from all firms that just defaulted, bearing on the banking sector. Note that 
< Po < P^- The parameter / € [0,1] reflects the impact of these defaults - either entirely on the cost of loans 
(/ = 1) or on the revenue of deposits (/ = 0). The logic behind rule Eq. (7) is that the interest rate on loan increases 
when defaults increase, while the rate on deposits is chosen in such a way that the profits of the banking sector are 
exactly zero at each time step. Indeed, one has, for any value of /: 

pHt)S-{t) - p'^{t)Xit) - V{t) = 0. ( 8 ) 

's ^-✓ V--✓ V—^ 

interest from loans interest paid on deposits cost of defaults 


Note that p’^{t) can become negative for large enough 'D{t) when / 7 ^ 1, i.e. the deposits are taxed by the banking 
sector. This could be interpreted as a kind of “bail-in tax” to absorb debt in extreme cases. In our simulations where 
/ = 0.5, this situation does occur in the unstable phases of the economy during a relatively short fraction of the time, 
corresponding to the peaks of the unemployment spikes. Finally, one could introduce an extra haircut in p^ or if 
one wants to model a profit-seeking banking sector, but the resulting profits would somehow have to be re-injected in 
the economy - an extra modeling step that we avoid at this stage by assuming the above no-profit rule. 


3. Households’ consumption budget 

As mentioned above, one major simplification of Mark-0 is to treat the whole household sector at the aggregate 
level, and to represent it with only a few variables: total savings S{t), total wages Wxit) = and total 

consumption budget Cs^t) (which, as emphasized in [I], is in general larger than the actual consumption C(t)). 

The effect of interest rates on households is two-fold. First, quite trivially, they receive some interest on their savings 
S{t) that adds to the wages Writ) and dividends as their total income. Second, the comparison between interest 
rates and inflation creates an incentive to consume or to save. This is in the spirit of the standard Euler equation of 
DSGE models where consumption is found to depend on the difference of rates on deposits p‘^{t) and inflation 7 r(t) 
(see e.g. [3, 35]). We therefore posit that the consumption budget of households C'B(t) is given by: 

Csit) = c{t) [S{t)+ WTit) + p‘^{t)S{t)] with c{t) = Co [1 + ad^t - pf)] , (9) 

where c{t) € [0,1] is the consumption propensity (taken to be a constant in Mark-0) and ac > 0 is a coupling constant 
that determines the sensitivity of households to the (moving average of the) difference between inflation and the 
interest paid on their savings. The larger the difference between the two, the larger the propensity to consume, as 
in standard equilibrium models [35], but with here undetermined phenomenological parameters co,ac that should in 
principle be measured on micro-data (surveys, laboratory experiments, etc., see e.g. [12-14]). 

Apart from these changes, the behavior of households is exactly the same as in Mark-0, see ]1] and Appendix A for 
details. 


4- Firms’ policy when debt is costly 

a. Financial fragility. Unlike households the N-p firms are heterogeneous and treated individually. Each firm is 
characterized by its production Yi (equal to its workforce), demand for its goods Di, price pi, wage Wi and cash 


® Note that the parameter / is similar to, but different from the parameter also called / in [1], which was used to share the cost of defaults 
on firms and households. 



balance Ei. The debt level of a firm is measured through the ratio 

= -E,l{W,Yi), 


( 10 ) 


which we interpret as an index of financial fragility. If $i(t) < 0, i.e. when the flux of credit needed from the bank 
is not too high compared to the size of the company (measured as the total payroll), the firm is allowed to continue 
its activity. If on the other hand $i(t) > 0, the firm defaults and contributes to total default costs T>{t). 

b. Production and wage update. If the firm is allowed to continue its business, it adapts its price, wages and 
productions according to reasonable “rules of thumb” introduced in [1] - see Appendix A. In particular, the production 
update is chosen as: 


liYi{t) < D,{t) => Yi{t + l)=Yi{t) + m\n{'qf{Di{t)-Yi{t)),u*{t)} 

Ifr,(t)>A(I) ^ Y,{t + l)=Y,{t)-p-[Y,(t)-D,{t)] 

where u*(t) is the maximum number of unemployed workers available to the firm i at time t (see Appendix A). The 
coefficients rj^ € [0,1] express the sensitivity of the firm’s target production to excess demand/supply, and they were 
constant in Mark-0. 

Here, we further postulate that the production adjustment depends on the financial fragility of the firm. Firms 
that are close to bankruptcy are arguably faster to fire and slower to hire, and vice-versa for healthy firms. In order 
to model this tendency, we posit that the coefficients rjf for firm i are given by: 

r]~ =? 7 omax(I-|-r$i(t), 0 ) 

? 7 += i?? 7 omax(I - r$i(t),0) (12) 

where rj^ is a fixed coefficient, identical for all firms, and R is the propensity ratio discussed in the previous section. 
The factor T > 0 measures how the financial fragility of firms influences their hiring/firing policy, since a larger value 
of then leads to a faster downward adjustment of the workforce when the firm is over-producing, and a slower 
(more cautious) upward adjustment when the firm is under-producing. The above definition however ensures that rj^ 
always remain non-negative, i.e. the reaction of the firms is always in the intuitive direction. 

It is plausible that the financial fragility of the firm also affects its wage policy: we give in Appendix A the wage 
update rules of Mark-0 and their modification to account for financial fragility, through the very same parameter F. 
In essence, deeply indebted firms seek to reduce wages more rapidly, whereas flourishing firms tend to increase wages 
quickly. 

The baseline Mark-0 model corresponds to F = 0, and leads to the phase diagram shown in Fig. I above. Interest¬ 
ingly, a non-zero value of F = Fq (constant across firms and in time) changes substantially the nature - but not the 
existence - of the phase transition between the full employment (FE) and full unemployment (FU) phase. The first 
order (discontinuous) transition for Fq = 0, 0 = oo found in [I] and shown in Fig. I, is replaced by a second order 
(continuous) transition when the firms adapt their behavior as a function of their financial fragility, i.e. when Fq > 0. 
Moreover, the “bad” FU phase for R < Rc becomes a partial unemployment phase with m < I that continuously 
varies with R: see Appendix B and Fig. 9 for full details. As firms become more careful, employment can be to some 
extent preserved - as expected. 

The “good” phase of the economy, on the other hand, is only mildly affected by a non zero Fq - for example the 
FE region of Fig. I expands downwards, which is expected since firms manage more carefully their balance sheet, 
reducing the occurrence of defaults. 

c. The influence of interest rates on the strategy of firms. We now argue that F should in fact depend on the 
difference between the interest rate and the inflation: high cost of credit makes firms particularly wary of going into 
debt and their sensitivity to their financial fragility should be increased. Therefore, we postulate that interest rates 
feedback into the behavior of the firm primarily through the F parameter, that we model in the simplest possible way 
as: 


F = max{ar(^(I) - 7r(t)),Fo}, (13) 

where or (similarly to olc above) captures the influence of the real interest rate f^if) — ^(t) on the hiring/firing policy 
of the firms. Whenever the real interest rate stays below Fo/orj the response of firms to changes of interest rates 
is negligible (perhaps as reported in [II]), whereas larger rates lead to a substantial change in the firms policy. The 
case Or = 0 but Fq > 0 corresponds to the above discussion and it is interesting in itself (see Appendix B). However, 
since we will be mostly concerned with policy issues, we will concentrate below on the other extreme case, ar > 0 
and Fq = 0, keeping in mind that reality is probably in-between. Note that F as defined above is now zero when real 
interest rates are negative, and is positive otherwise. 
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C. Summary 

1. How many new assumptions? 

In our attempt to include interest rates and monetary policy into the Mark-0 framework, we have made new 
behavioral assumptions, that we tried to keep as simple and as parsimonious as possible: 

• For the Central Bank, we have merely adapted the standard Taylor rule to our setting, introducing three 
monetary policy parameters p* (the natural interest rate) and 4>Tr,4>e (the Taylor rule parameters), and two 
targets: inflation tt* and employment e*. 

• For the banking sector, we have essentially used trivial accounting rules to ensure a no-profit condition on loans 
and deposits. The only new parameter is / and determines how much of the cost of bankruptcy must be paid by 
loans or by deposits. However, this parameter plays very little role in the qualitative behaviour of the economy, 
so we set it to /= 1/2 henceforth. 

• For the households, we assume an Euler-like behaviour of the consumption budget as a function of the real 
interest rate on deposits, i.e. the consumption budget decreases when the real rate is high and increases when 
it is low, with a slope given by parameter Oc- 

• For the Arms, the hiring/firing and wage policy is affected by the real rate on loans. When it is high, indebted 
firms are more careful, firing more rapidly and hiring less. The unique parameter coupling the real interest rate 
to the firms policy is or- 

These last two additional rules are to some extent arbitrary. However, they capture effects that certainly exist in the 
real world; furthermore changing the detailed implementation we chose here while keeping the spirit of these rules 
lead to very similar conclusions. This is in line with our general claim that macro-behaviour is to a large degree 
insensitive to micro-rules. 

In summary, we want to emphasize that our extension of Mark-0 is very parsimonious: we have only added two 
behavioral parameters {ac,ar)- 


2. Recovering Mark-0 

From the above discussion, we see that the core Mark-0 model of [1] is recovered whenever the baseline interest 
rate is zero p* = 0, the CB is inactive ((/,r = i's = 0), and households and firms are insensitive to interest rates and 
inflation (i.e. setting «(, = c^r = 0). 

There is however a slight remaining difference with Mark-0 in the resolution of bankruptcies. The closest one can 
get is by setting / = 0, i.e. absorbing default costs only through savings. In this case the only non-zero interest rate 
remaining in the dynamics of the model is the one on deposits, which is negative: p'^ = —'D{t)/X{t) < 0. This indeed 
roughly corresponds to the default resolution described in [1] where default costs are paid by households and firms 
savings.® There are also minor differences in the time-line of the model (in particular bankruptcies are resolved before 
price, production and wages are updated). All these differences however have a negligible quantitative impact on the 
results below. 


III. THE “NATURAL” BEHAVIOR OF THE ECONOMY (WITHOUT MONETARY POLICY) 

In this section we analyze the features of the model that arise from the introduction of interest rates in the Mark-0 
economy, disregarding for a while any active monetary policy (i.e. setting /).„. = = 0 in Eq. (3) above). In other 

words, we study an economy where the baseline interest rate is equal to p*, constant in time, and affects both the 
consumption propensity of households through the parameter ac appearing in Eq. (9), and the firms hiring/firing 
propensity through the parameter op appearing in Eq. (13), with Fq = 0 henceforth. 

When interest rates do not feedback at all into firms’ and agents’ behavior (i.e. for ap = Tq = 0 and Oc = 0) the 
phenomenology of the Mark-0 model is basically unaffected; in particular the phase diagram Fig. 1 is unchanged. 


® To be more precise, in the default resolution described in [1] we introduce a bailout probability, called / there, which sets the relative 
impact of default costs on households and firms savings. In this sense, the present setting recovers the one in [1] with / 1/2. 
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FIG. 2: Left: Phase diagram of the “natural” state of the economy in the p* — 0 plane (i.e. with no monetary policy, 
(j>^ = = 0) when both firms’ and households’ decisions are sensitive to the rate level: ap = 50 and = 4. For small enough 

rates p* < p** , one recovers the FE-EC-RU phases as 0 is decreased. When p* > p** , however, the RU phase prevails for all 
values of 0, which becomes irrelevant. Right: Phase diagram in the p* — ar plane for 0 = 3, showing the dependence of the 
critical value p** on ap. The other parameters of the model are set to: R = 2 (with 77 ° = 0.1), co = 0.5, Fq = 0, v? = 0.1, 
7 p = 7 m = 0.05, /3 = 2, 5 = 0.02, / = 0.5, Nf = 2000 (see Appendix A for a definition of all these parameters). 


A. Coupling between interest rates and firms behavior 

When ar > 0, on the other hand, the overall behavior of the economy evolves as expected: as long as p* is less than 
the inflation rate, nothing much happens, in particular because Eq. (13) gives E = Eq = 0 here. When p* exceeds the 
inflation rate, one observes that the unemployment rate u starts increasing with p*, while the demand for credit and 
the inflation rate itself nosedive as expected. 

In Fig. 2 (left) we plot the phase diagram of the model in the p* — 0 plane when R = 2 > Rc, and for a fixed value 
ar = 50.^ For p* smaller than a certain value p** « 1.3, one observes the familiar three phases FE-EC-RU as 0 is 
decreased, as in Fig. 1. However, for a baseline rate larger than p**, the FE and EC phase disappear entirely, and 
the Residual Employment phase (with u ^ 30%) prevails for all values of 0. To wit, when interest rates are too high, 
firms hesitate to accumulate more debt (even if they are allowed to when 0 is high). Rather, they prefer reducing 
their work force when needed, keeping the unemployement at a relatively high level. 

Fig. 2 (right) allows one to understand the role of ap: there we show the phase diagram in the p* — ap plane for 
fixed values of i? = 2 and 0 = 3, such that the economy is in the Full Employment phase for ap = 0. A sudden phase 
transition between the FE and RU phases occurs for a critical value p**(ap); the larger the sensitivity to interest 
rates - i.e. the larger ap - the smaller the critical value of the baseline interest rate beyond which the economy is 
destabilized. It is interesting (and quite counterintuitive) that the aggregate behavior of the economy is not a smooth 
function of the interest rate. When Arms are risk averse and fear going into debt, large enough interest rates lead to 
more unemployment that spirals into a destabilizing feedback loop. This is one of the “dark corners” that ABMs can 
help uncovering. 


B. Coupling between interest rates and household behavior 

Perhaps unexpectedly, the coupling between interest rate and consumption (captured by parameter ac) appears to 
have much smaller influence on the “natural” state of the economy, at least when ac is chosen within a reasonable 
range. Its main influence is to increase the output fluctuations around the steady state, by amplifying price trends 
through the resulting reduction/increase in consumption. Interestingly, we find that for 0^1 and independently of 


^ This value of op has the following interpretation: when the debt of a firm equals its payroll, i.e. when = 1, a real interest rate of 8% 
annual leads to a freezing of all hires and a doubling of the firing rate, compared to a zero-debt situation. 
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FIG. 3: Aggregate response of the economy to a policy shock represented by a drop in the interest rate po set by the central 
bank (from 2% to 1.8%). The vertical dashed line gives the timing of the policy shock. All aggregate variables are represented 
as a relative variation from their average values prior to the shock. For a better visualization, we amplified by a factor 100 the 
wage and price responses. Note that the damped oscillation relaxation towards equilibrium. Parameters are 0 = 3, p* = 2%, 
ar = 50, Qc = 4, and other parameters as in Fig. 2. 


p*, ac has a stabilising effect on the economy: i?c shifts to lower values as Oc increases (in the absence of monetary 
policy). Clearly, micro-data is needed to estimate the value of Oc for realistic applications; but since this parameter 
plays a small role we will choose rather arbitrarily ac = 4, unless explicitly stated. This corresponds to a moderate 
sensitivity to inflation/interest rates: a rise of the interest rate of 1%/year corresponds to a decrease of 4% of the 
consumption propensity. 


IV. MONETARY POLICY EXPERIMENTS 

We now consider a simple policy framework where the Central Bank adjusts the base interest rate po(t) in order 
to achieve its inflation and employment targets. Given the simplicity of our model we are mainly interested here in 
gaining a qualitative understanding of the possible consequence of a Taylor-rule based monetary policy in a stylized 
Agent-Based framework. We defer parameter calibration, detailed comparison with simple DSGE models and more 
quantitative insights to future studies. 

Our major finding is that provided its policy is not too aggressive and the economy not too close to a phase 
transition, the CB is successful in steering the economy towards its targets. However, the mere presence of different 
equilibrium states of the economy separated by phase boundaries (i.e. “dark corners”) may deeply alter the impact 
of monetary policy. Indeed, we will exhibit cases where the monetary policy by itself triggers large instabilities and is 
counter-productive. 

As in the previous section, we will refer to the state obtained without any response of the CB, i.e. when 4>Tr = 4>e = Oj 
as the “natural” state of the economy (for a given set of parameters). The corresponding “natural” value of a variable 
X will be denoted by Xnat- In order to simplify the analysis and since most of the parameters of our model play little 
role in the qualitative behavior of the economy we set once and for all some of them to the values given in the caption 
of Fig. 2 and choose p* = 2%. We only focus on the four parameters defining the CB policy (i.e. (fej £* and 
TT*), the parameters of the transmission channels (i.e. ac and ar), and the two parameters locating the system in the 
phase diagram of Fig. 1 (i.e. the hiring/firing ratio R and the bankruptcy threshold 0). 


A. A monetary shock 


In order to check that our framework produces reasonable results, it is interesting to start with the simplest policy 
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FIG. 4: In these plots we select the same natural state of the economy with Q = 2, R = 2, p* — 2% and ar = 50, such that 
rtnat ~ 0.33 and TTnat ~ 0 (but fluctuating). The CB sets a much lower unemployment target of u* = 0.05 and an inflation 
target of tt* = 0.2%. In the left plot the CB policy is moderate = cj>£ = 0.5) and basically achieves its goals, while in the 
right plot the policy is aggressive (<)>,r = <)>£ = 1) and destabilizes the whole economy. In the insets, po and tt are given in %s. 
The period of the business cycles in the right plot corresponds to « 4 years if the elementary time unit of the model is one 
month. Larger spikes of unemployment are also observed with much lower frequency (not shown). 


experiment, i.e. an expansionary monetary policy shock where the baseline interest rate is instantaneously decreased 
from 2% to 1.8%, all other parameters of the model being fixed, and with no further intervention of the Central Bank 
(i.e. (^,r = = 0). The resulting dynamics of output, wages and prices is shown in Fig. 3 and can be compared 

with e.g. Fig 1 of Ref. [36]. It is gratifying to see that our simple modelling strategy, in particular the way the 
monetary policy is channelled to the economy, does lead to a quite realistic shape of impulse response functions, at 
least qualitatively. In fact, the empirical hump-shaped response of the output reported in [36] is quite well reproduced 
by their DSGE model with nominal rigidities. Although the underlying microfundations are very different, our ABM 
indeed includes frictions: all adjustments (of production, wages and prices) in our model are only made progressively. 
So our framework is indeed able to capture effects predicted by enhanced DSGE models. No attempt was made here 
to adjust parameters to quantitatively match empirical data, in particular the time scales: this is beyond the scope 
of the present study which is primarily intended to be a description of the qualitative behavior of our toy-economy. 


B. Mild vs aggressive monetary policy 

In Fig. 4 for example, we show the result of the policy of the Gentral Bank that attempts to bring down the natural 
unemployment rate of Unat ~ 0.33 (a rather large value corresponding to 0 = 2, i? = 2 and p* = 2%) to a low target 
of I — e* = M* = 0.05. The target inflation is tt* = 0.2% per time step (corresponding to 2.4% annual if one interprets 
the time step to be a month), compared to a natural inflation that fluctuates around zero: Tr^at ~ 0. The left graph 
corresponds to a mild monetary policy, with Taylor-rule parameters set to = (p^ = 0.5. The policy is seen to be 
rather successful: the inflation is on target, while unemployment goes down to m « 0.07, not far from the target of 
0.05. But now look at the graph on the right, where the only difference is the aggressiveness of the GB that attempts 
to reach target too quickly, merely doubling the value of I. In this case, the monetary policy has induced 

strong instabilities, with “business cycles” of large amplitude and inflation all over the place. 

We show in Fig. 5 the result of an extensive exploration of the role of and pg when 0 = 2, i? = 2 and 
p* = 2%. One sees that there is a wedge-like region around = 0 where the policy does not induce instabilities 

(signaled by a yellow/orange hue in Fig. 5-right). However, the region of parameters where the unemployment rate is 
significantly reduced is only a subset of this wedge, corresponding to the black region in Fig. 5-left, where pg ~ 0.5, 
(pT^ < 0.5. In other words, the Gentral Bank must navigate in a narrow window: too little is not enough, too aggressive 
is counterproductive and leads to instabilities and wildly oscillating economies. 

Fig. 5 also reveal that in our toy-world, a Gentral Bank with a dual mandate (inflation and output) enjoys a wider 
region of stability compared to a Gentral bank with inflation as the only mandate. Steering inflation and output 
simultaneously increases the probability of a successful monetary policy. 
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FIG. 5: Policy performance in the {(pn, (pe) plane (the origin corresponds to the natural state of the economy), for O = 2, R — 2, 
p* = 2%, ar = 50 and Qc = 4. The target inflation and unemployment are, respectively, tt* = 0.2% and u* = 0.05. Left: Color 
map of the average unemployment; Right: Color map of the amplitude of the business cycle, measured as maxt(M) — mint(u). 
Yellow/orange regions correspond to unstable economies with crises of large amplitude. As one can see the policy is effective 
as long as it is not too aggressive, with a sharp transition to a regime where it becomes detrimental. The closer the natural 
economy to a phase boundary, the more destabilizing the policy - see Fig. 6. 


C. The role of phase boundaries 

The fragility of the economy is clearly due to the proximity of the phase boundary that appears in Fig. 2. As one 
moves away from the boundary, for example by increasing 0, one finds that the region where the CB policy is harmful 
shrinks. We illustrate this by moving along the line cpT^ = (pe = 4>cb in parameter space. We display in Fig. 6 the 
phase diagram of the model in the (0 cb, 0) plane and in the {(PcbtR) plane. One clearly sees from the left graph 
on the top row that the deep blue region (corresponding to low unemployment) expands for larger 0, and that the 
yellow/orange region of the graph on the right (corresponding to strong oscillations) recedes. 

The bottom graphs illustrate the role of R. One mostly observes that: 

• (a) deep in the FU phase (i? < 0.75), the monetary policy is helpless in restoring employment; 

• (b) for intermediate values of R, large enough values of ^cb do lead to small unemployment rates; 

• (c) when R and (pcB are simultaneously large, instabilities appear (see Fig. 6 right bottom graph, north-east 
corner). 

Point (a) above is interesting and can be understood as follows: when R is small, firms are so quick to adjust production 
downwards that they never need credit, and their financial fragility is low or even negative. The interest-rate impact 
parameter ar then becomes completely ineffective in this case. This could be relevant to understand the aftermath 
of the 2008 crisis: if one interprets strong downsizing (i.e. small R) as a result of a drop of confidence induced by the 
Lehman crisis, the above discussion suggests that a low interest rate policy might not be as effective as one may have 
expected. 

As a complement to the above analysis we also investigated the performance of the policy as a function of the 
distance between target and natural values. We find, not surprisingly, that when targets are not too far from the 
natural state of the economy the CB manages to achieve its targets without triggering any instability. When however 
targets are far from the natural state even a mild policy may become detrimental and trigger instabilities. 


D. The role of households and firms sensitivity to rates 

We have also investigated the role of the two policy transmission channels (ar, firms and Oc, households) separately, 
and found that both channels in isolation may trigger instabilities. In Fig. 7 we plot the policy performance in the 
(4>CB, cuc) plane and the ((()CB, cir) plane, with all other parameters fixed, in particular R = 2 and 0 = 3. As expected, 
the top row shows that the larger the value of ac, the more careful the monetary policy has to be in order to avoid 
instabilities. Note, interestingly, that there is a thin region in the {(pCB,o.c) plane (spot the yellow dots) where 
unemployment goes to u « 1, i.e. the economy is completely destabilized by the monetary policy! 







14 



FIG. 6: A more systematic diagram of the effect observed in Fig. 5. Here we move on the diagonal <j)-n = (f)s of Fig. 5 on the 
(?!>CB-axis and plot the average unemployment u and unemployment variations maxt(M) — mint(«). The other axis is either 0 
(top row, for R = 2) or R (bottom row, for 0 = 10), still with ar = 50, ac = 4., tt* = 0.2% and u* = 0.05. Note that since 
Or > 0, the transition at Rc « 1,<(>cb = 0 is second order (see Appendix B). As above, yellow/orange regions in the plots on 
the right, correspond to unstable economies with crises of large amplitude. 


The bottom row of Fig. 7 shows the interplay between policy aggressiveness c^cb and firms sensitivity to the real 
interest rate ap- When c^cb = 0, one recovers the FE-RU transition for ap ~ 30 already observed in Fig. 2. For 
larger op’s, the CB policy is at first successful in reinstalling the FE phase, before destabilizing again the economy 
beyond « 0.9. When op becomes small, a portion of the FE region expands up to higher values of ^cb ~ 1-25. 

A detailed comparison between Fig. 6 (top row) and Fig. 7 (bottom row) reveals that the simultaneous presence of 
the two transmission channels has an overall stabilizing effect. Indeed, when only one channel is present (i.e. either 
Oc or Op is zero) the value of (/>cb beyond which the system is unstable is smaller than when both channels are active. 
For example, for 0 = 3, Oc = 4. and op = 50. one sees from Fig. 6 (top row; left, horizontal line) that the economy is 
destabilized for (j)^^ « 1.3, to be compared to « 0.9 when ac = 0 and to 4>q-q « 0.5 when ap = 0. Although we 
cannot provide a complete interpretation, this interesting observation might be due to the fact that the a parameters 
can be interpreted as “awareness” parameters that allow firms and agents to better anticipate crises and thus dampen 
the destabilising influence of the Central Bank policy. 


E. Comparison with DSGE models 

It is quite interesting to compare the diagram of Fig. 5 with its DSGE counterpart (see e.g [3], section 4.3). Within 
our Agent-Based framework, small values of are at worst ineffective wheras large values of (j)-K,4‘e lead to 

instabilities. This is at the opposite of what is observed in the fully rational world of DSGE models, where small 
values of 0^, lead to instabilities, while large values of allow the Central Bank to funnel the economy on a 

stability path [3]. The intuition for this is made crystal clear by Gali in [3]: The monetary authority should respond 
to deviations of inflation and the output gap from their target levels by adjusting the nominal rate with “sufficient 
strength”; [...] it is the presence of a “threat” of a strong response by the monetary authority to an eventual deviation 
of the output gap and inflation from target that suffices to rule out any such deviation in equilibrium. In other words. 
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FIG. 7: Policy performance, as in the previous two plots, now separately as a function of the households sensitivity to rates 
ac and the firms sensitivity to rates ar- The i-axis is again = 4>cb., whereas Q = 3, R — 2 and p* = 2%, with 

again tt* = 0.2% and u* — 0.05. (Top row) Here ar = 0 and ac > 0. The FE (black) region prevailing for = (j>cB = 0 is 
destabilized beyond a certain <(>qb that decreases as increases. Note the constellation of yellow points that correspond to 
a completely broken economy. {Bottom row) We now set = 0 and vary op > 0. When (j>cB = 0, one recovers the FE-RU 
transition for ar ~ 30 observed in Fig. 2. One sees that for larger ar’s, the CB policy is at first successful in reinstalling the 
FE phase, before destabilizing again the economy beyond the vertical line <(>qb « 0.9. 


fully rational, forward looking agents know that inflation will be tamed by the response of central authorities, an 
element that is indeed completely absent in the myopic forecast world of most ABMs (including Mark-0), where 
agents only use their knowledge of the past and present situation and adapt their behavior accordingly. The fact that 
the outcome of these two hypotheses on the stability of the economy are so radically different should be a strong caveat. 
Which of these two extreme visions of the world is closest to reality is of course a matter of empirical investigation 
(on this point, see [37], and for recent attempts to include deviations from perfect rationality in DSGE models, see 
e.g. [38]). 


V. SUMMARY, CONCLUSION 

In this paper, we parsimoniously extended the stylized macroeconomic Agent-Based Mark-0 model introduced 
in [1], with the aim of investigating the role and efhcacy of monetary policy. We focused on three effects induced by a 
non-zero interest rate in the model, that are arguably the most important transmission channels of the Central Bank 
policy: i) a change of the accounting rules to factor in the cost of debt and the extra revenue of deposits; ii) a change 
of the consumption behavior of household that depends negatively on the real interest rate and iii) a change in the 
hiring/firing and wage policies of firms, that avoid running into debt when interest rates increase. This amounts to 
adding only two new parameters to the baseline Mark-0 model of [1]. 

We first studied the model in the absence of monetary policy, i.e. without a “Taylor-rule” that creates a feedback 
between inflation, unemployment and interest rates. The introduction of a coupling between the financial fragility of 
firms and the hiring/firing and wage policies has two main effects: a) the first order (discontinuous) phase transition 
between a ‘good’ and a ‘bad’ phase of the economy, discussed in [1] is replaced by a second order (continuous) 
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transition; b) a new first order transition to a phase with high residual unemployment (RU) appears when the interest 
rate is larger than some critical value, even in the region where full employment (FE) is achieved for zero-interest rate 
(i.e. 0 ^ 1 in Fig. 1). The larger the sensitivity to interest rates, the smaller the value of the baseline rate beyond 
which the economy is destabilized. In that region, allowing firms to accumulate more debt does not help stabilizing 
the economy. 

We then allowed the Central Bank to adjust the baseline rate so as to steer the economy towards prescribed levels 
of inflation and employment. Our major finding is that provided its policy is not too aggressive (i.e. when the 
targets are not too far from the ‘natural’ state of the economy, and for a low enough adjustment speed) the Central 
Bank is successful in achieving its goals. However, the mere presence of different states of the economy separated by 
phase boundaries, besides being interesting per se, can cause the monetary policy itself to trigger instabilities and be 
counter-productive. The destabilizing influence of the Central Bank also depends on the firms/households sensitivities 
to rates. Perhaps ironically, too small sensitivities make the Central Bank policy inefficient, but too large sensitivities 
make the same policy dangerous. Seen differently, the Central Bank must navigate in a narrow window: too little is 
not enough, too much leads to instabilities and wildly oscillating economies [39]. As mentioned in the last paragraph, 
this conclusion strongly contrasts with the prediction of DSCE models. Interestingly, we also find that a Central 
Bank with a dual mandate (inflation and output) enjoys a wider region of stability compared to a Central bank with 
inflation as the only mandate. 

As we emphasized in the introduction, the key message of both our previous paper [1] and the present one is that 
even over-simplified macroeconomic ABMs generically display a rich phenomenology with an economy characterized 
by different states separated by phase boundaries across which radical changes of the emergent behavior take place. 
These are, we believe, the “dark corners” alluded to by O. Blanchard in [5], that both academics and policy makers 
should account for and wrestle with. We believe that the major advantage of ABMs over DSGE-like models is the 
very possibility of crises at the aggregate level, mediated by generic feedback mechanisms whose destabilizing role 
may not be immediately obvious or intuitive. Rather than the precisely calibrated predictive tools that standard 
equilibrium models claim to provide [24], ABMs offer extremely valuable qualitative tools for generating scenarios, 
that can be used to foresee the unintended consequences of some policy decisions. Some of this outcomes, which would 
be “Black Swans” [40] in a DSGE framework, can in fact be fully anticipated by schematic ABMs [41]. As expressed 
with remarkable insight by Mark Buchanan [42]: Done properly, computer simulation represents a kind of “telescope 
for the mind, ” multiplying human powers of analysis and insight just as a telescope does our powers of vision. With 
simulations, we can discover relationships that the unaided human mind, or even the human mind aided with the best 
mathematical analysis, would never grasp. 
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Appendix A: Price, production and wage updates 
1. The households timeline 

• At the beginning of the time step households are characterized by a certain amount of savings S(t) > 0. 

• After each firm chooses its production Yi{t), price Pi{t) and wage Wi{t) for the current time step (see later) 
wages are paid. Since firms use a one-to-one linear technology taking only labor as input and productivity is 
set to 1, the production equal the workforce of the firm. Hence, the total amount of wages paid is given by 

WTit) = J2w^{t)Y,it) (Al) 

I 

• Once the total payroll of the economy is determined, interests on deposits are paid and households set a 
consumption budget as 


CB{t) = c{t)[S{t) + Writ) + /(t)^(t)] 


(A2) 



where c{t) G [0,1] is the propensity to consume and may depend on inflation/interests on deposits, see Eq. (9). 

The consumption budget is distributed among firms using an intensity of choice model [43]. The demand of 
goods for firm i is therefore: 

where j3 is the price sensitivity parameter determining an exponential dependence of households demand to the 
price offered by the firm. Indeed, /3 = 0 corresponds to complete price insensitivity and /3 —)■ oo means that 
households select only the firm with the lowest price. In this sense, as long as /3 > 0 firms compete on prices. 

The actual consumption C(t) (limited by production) is given by 

Afp Nf 

C{t) := p^{t) min {Y,{t),Di{t)} < CbH) = (A4) 

and households accounting therefore reads 

Sit + l) = [l + /(t)]5(t) + IET(t)-^a(t)min{A,Ej + A(<) (A5) 

i 

where A{t) are dividends paid (see below for a definition of this last term). 


2. The firms timeline 


At the beginning of the time step t firms with 4>i(f) > 0 become inactive and are removed from the simulation. 
Default costs are computed as 

m = - E (A6) 

i bankrupt 

Firms with 4>i(t) < 0 are instead allowed to continue their activity and contribute to total loans £~{t) and 
total firms savings £~^{t) as 

£- = - ^ min{£i(f),0} 

i not bankrupt 

= E max{£li(t),0}. (A7) 

i not bankrupt 


Active firms set production, price and wage for the current time step following simple adaptive rules which are 
meant to represent an heuristic adjustment. In particular: 


- Price: 

Prices are updated through a random multiplicative process which takes into account the production- 
demand gap experienced in the previous time step and if the price offered is competitive (with respect to 
the average price). The update rule for prices reads: 


If y,(f) < A(t) 


ify,(f)> A(t) 


<p{t) 
lipt{t) >p{t) 

liPt{t) >p{t) 
lipt{t) <p{t) 


Pi{t+l) =p,(f)(I -k7p^,(t)) 

Pi{t -b 1) = Pi{t) 

Pi{t + l) =p,(t)(I - 7p^,(t)) 

Pi{t -b 1) = Pi{t) 


(AS) 


where ^i{t) are independent uniform {7[0,1] random variables and 7 p is a parameter setting the relative 
magnitude of the price adjustment (we set it to 5% unless stated otherwise). Fig. 8, which plots the 
average profit of firms as a function of the offered price, shows that these rules lead to a reasonable 
emergent “optimizing” behavior of firms. As expected, the profit reaches a maximum for prices slightly 
above the average price p{t) . Higher prices are not competitive and firms lose clients, while lower prices do 
not cover production costs. 
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FIG. 8 : Scatter plot in one time step of firm profits vs the price offered. The red line correspond to a moving flat average of 
100 consecutive points. Parameters are: i? = 2 (with 77 ° =0.1), 0 = 2, 7 ^ = 7 „ = 0.05, ar = 50, Fq = 0, /3 = 2, = 4, 

N = 5000 and p* = 0.5%. 


- Production: 

Independently of their price level, firms try to adjust their production to the observed demand. When 
firms want to hire, they open positions on the job market; we assume that the total number of unemployed 
workers, which is N-pult), is distributed among firms according to an intensity of choice model which 
depends on both the wage offered by the firm® and on the same parameter /3 as it is for Eq. (A3); therefore 
the maximum number of available workers to each firm is: 

^PWi{t)/w{t) 

^ ^ ^l3Wi(t)lw(t) -^Fu(t) ■ (A9) 

The production update is then defined as: 

< Diit) ^ l)(t + l) =l)(t)+min{?7+(A(t)-h)(t)),M*(t)} 

liY,{t)> D,{t) ^ Y,{t+l)=Y{t)-T^-[Y,{t)-D,{t)] 

where 77 ^ S [0,1] are what we denote as the hiring/firing propensity of the firms. According to this 
mechanism, the change in output responds to excess demand (there is an increase in output if excess 
demand is positive, a decrease in output if excess demand is negative, i.e. if there is excess supply). The 
propensities to hire/fire r]± are the sensitivity of the output change to excess demand/supply. 

- Wage: 

The wage update rule we chose follows (in spirit) the choices made for price and production. We propose 
that at each time step firm i updates its wage as: 


wlit +1) = w,(t)[i + 7v,(i - r$,)£?'(t)] 

w,{t +1) = w,(t)[i - 7»(1 + rd>o<W] 

where 7 ^, is a certain parameter, Vi{t) is the profit of the firm at time t and ^'(t) an independent f7[0,1] 
random variable. If {t + 1) is such that the profit of firm i at time t with this amount of wages would 
have been negative, Wi{t + 1) is chosen to be exactly at the equilibrium point where Vi(t) = 0 otherwise 
W,{t + l) = Wf{t + l). 

The above rules are intuitive: if a firm makes a profit and it has a large demand for its good, it will 
increase the pay of its workers. The pay rise is expected to be large if the firm is financially healthy and/or 


if 


if 


n{t) 

Ym 

V.(t) 


< Ui{t) 

> 0 

> D,{t) 

< 0 


(All) 


A higher wage translates in the availability of a larger share of unemployed workers in the hiring process. 
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if unemployment is low (i.e. if e is large) because pressure on salaries is high. Conversely, if the firm makes 
a loss and has a low demand for its good, it will reduce the wages. This reduction is drastic if the company 
is close to bankruptcy, and/or if unemployment is high, because the pressure on salaries is then low. In all 
other cases, wages are not updated. 

The parameters allow us to simulate different price/wage update timescales, i.e. the aggressivity with 
which firms react a change of their economic conditions. In the following we set 7 ^ = 0.05 and = 7 ^. 
The case 7 iu = 0 corresponds to removing completely the wage update rule, such that the version of the 
model with constant wage is recovered. 

• After prices, productions and wages are set and interests paid, consumption and accounting take place. Since 
each firm has total sales Pi min {!/, firms profits are 

'Pi{t) =pj(t)min{i;(t), A(t)} - Wi{t)Y^{t) + max{£:i(t), 0} + min 0} . (AI2) 

When firms have both positive Si and Vi dividends are paid as a fraction <5 of the firm cash balance Si. The 
update rule for firms cash balance is therefore 

S,{t + 1) = S,{t) + V^{t) - 5S^{t)6{V^{t))6{S,{t)) (AI3) 

where 9(x) = 1 if x > 0 and 0 otherwise. Correspondingly, households savings are updated as 

S{t + 1) = S{t) + ^ - '^pi{t) min {Yi{t),Di{t)} + 6'^S^{t)0{Vi{t))0{Si{t)). (AI4) 

i i i 

The dividends share <5 is set to 2% unless stated otherwise and the A(t) term in Eq. (A5) is given by 

A(t) = bY,S,{t)0{VmK^^^)) (AI5) 


• Finally, an inactive firm has a finite probability p (which we set to 0.1) per unit time to get revived; when 
it does so its price is fixed to Pi{t) = p{t), its wage to Wi{t) = w{t), its workforce is the available workforce 
Yi{t) = pu{t) and its cash-balance is the amount needed to pay the wage bill Si{t) = Wi{t)Yi{t). This small 
’liquidity’ is provided by firms with positive Si in shares proportional to their wealth Si. 


Appendix B: Firms’ adaptive behavior leads to a second order phase transition 


We start by analysing the effect of adaptation of firms. In order to get a first insight it is useful to consider a 
simplified setting where T = Tq (i.e. ar = 0), p^(t) = Po(i) = 0, / = 1, c{t) = Cq = 0.5 (hence Oc = 0) and wages are 
constant and equal to 1 {"fw =0). In this case the basic model described in [1] (with constant wages) is recovered. 

Intuitively, the coupling between financial fragility and hiring/firing propensity should have a stabilizing effect on 
the economy. Moreover, the full unemployment phase at R < Rc is deeply affected by the presence of ar: for ar 7 ^ 0 
the unemployment rate in this phase is no longer one, but becomes smaller than one and continuously changing with 
R. In order to derive an estimate of these continuous values we use an intuitive argument (at 0 = 00 ) which is 
justified a posteriori by the good match with numerical results. Given the critical ratio R = = Rc separating 

the high/low unemployment phases when there is no adaptation (i.e. Tq = 0 ) one can expect that equilibrium values 
of the unemployment rate different from 0 and 1 can only be stable if rfj^/rf_ remains around the critical value Rc at 
To = 0. Near criticality therefore we enforce that: 


^ ^ ??+(!-ro4>,) ^ _r * ~ - 0 + 

riS ?7°(I-hro$i) ° RcV-+'n+' 


(Bl) 


An explicit form of $ in terms of the employment rate e = Y can be obtained with the additional assumption that 
the system is always close to equilibrium (i.e. p ~ 1 and U ss F, at least when r]'\_/rf_ ~ i?c), which allows one 
to express households savings in terms of the firms’ production. Indeed (see the discussion in [I]) at equilibrium 
W = B = N-pY = c{W + S), from which it follows that NpY = W = >S'c/(I — c). For c = 0.5 as in our simulations one 
thus has S = NpY. Since the total amount of money is conserved (in our simulations NpS + S = NpS + NpY = Np, 
see Appendix C) one finally obtains that L = \ — Y and $ = (F — I)/F = (e — !)/£:, hence 

Fq _ Ret- -0% 

e ReP- +V+ ° Rc + R ° 


(B2) 
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FIG. 9: Inverse of the average employment rate To/e as a function of the ratio R = ri%/rf_ with rf/_ — 0.1 and 7 = 0.01 when 
To > 0. When the employment rate is rescaled with the parameter Fq (here Fq = 10“®, 10“^) the different lines collapse 
and Eq. (B2) agrees with numerical simulations. In the inset we also plot the rescaled variance, still as a function of 
Approaching the critical point the variance of the unemployment fluctuations diverges, together with their relax;ation time 
going to inhnity. The other parameters are: 5 = 0.02, 0 = 5, 7 ™ = 0, c = 0.5, /3 = 0 and ip = 0.1 


Note that according to this formula the employment goes to £ = 1 at the critical point R = Re- Above this value, 
the economy is in the “good” state and employment sticks to e = 1 (this is because in the argument the effect of 0 
has been neglected). Moreover, when R = Rc, £ is proportional to Fq and therefore in the limit Fq —>■ 0 one has £ = 0 
for all R < Re- This is the “bad” phase of full unemployment at Fq = 0, which becomes in this case a phase where 
employment grows steadily but remains of order Fq except very close to the critical point. 

Eq. (B2) is plotted in Fig. 9 together with numerical results. Note that in this case the representative firm 
approximation (iYp = 1) is in good agreement with numerical results also for Np = 10,000, as it was for the 
discontinuous transition obtained for F = 0. In the inset of Fig. 9 one can see that the variance of the fluctuations of 
employment rate is diverging as long as the critical value of R is approached. This is confirmed by a spectral analysis 
of the unemployment time series (see Fig. 10). In order to obtain the power spectrum we apply the GSL Fast Fourier 
Transform algorithm to the time series £{t) — (e). As one can see in Fig. 10 the power spectrum is well approximated 
by an Ornstein-Uhlenbeck form: 


I{uj) — Iq 


LOn 


■ OJ^ 


(B3) 


with Wo going linearly to 0 when rj^ approaches its critical value, meaning that the relaxation time diverges as 
one approaches the critical point. Note that this is not the case for the Mark 0 model with Fq = 0 which instead has 
a white noise power spectrum even in proximity of the transition line. The first order (discontinuous) transition for 
Fo = O,0 = oois thus replaced by a second order (continuous) transition when the firms adapt their behavior as a 
function of their financial fragility. 

Finally, note that the presence of a continuum of states for the unemployment rate whenever Fq > 0 and R < Rc 
holds also with 7 uj > 0 (when wages are not constant). It was however simpler to perform analytical computations 
with constant wages. 


Appendix C: Pseudo-code of Mark 0 

We present here the pseudo-code for the Mark 0 code described in Sec. IIB and Appendix A. The source code of 
the baseline Mark-0 is available on demand. 
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FIG. 10: Logarithm of the normalized power spectrum for Mark 0 with adaptive firms (Fq = 10“^), 7 p = 0.05 and A^f = 1 
(left) and Np = 1000 (right). The other parameters are set as in Fig. 9. The main plot shows two examples of the spectrum for 
jyK. =0.1 and 77 ° = 0.05 (black line) and 77 ° = 0.09 (red line) in the left plot, 77 ° = 0.05 (black line) and = 0.08 (red line) in 
the right plot. The time series is made of 2^® time steps after T^q = 500 000 and the logarithm of the spectrum is averaged over 
a moving window of 100 points for a better visualization. With both system sizes the fit with Eq. (B3) (blue dashed lines) is 
good with the only difference that when A^p > 1 a clear oscillatory pattern appears at high frequencies, that becomes sharper 
and sharper as A^f increases. In the inset of each hgure we plot the value of ujq in Eq. (B3) obtained from the ht as a function 
of the ratio R — rj^/rj^. In both cases cuo goes linearly to 0 as the critical value is approached. 
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Algorithm 1 Mark 0 


Require: Nf Number of firms; co,/i, 7p,7m, 


ri%,ri^, 5, 0, if, /, Oc, Qr, To, £*, tt*, p* ,ijJ\ T total evolution time; 

0 Initialization 


for (i 0; i < Nf \ i <— i + 1) do 
W[{] ^ 1 

p[i] ■<— 1 + 0.2(random — 0.5) 

Y[i] [1 + 0.2(raiidom — 0.5)]/2 
D[{\ ^ 0.5 

S[i] t— kF[i]y[z] random 

V[i] ■<— p[i] min(D[i], y[i]) — iy[i]y[i] 


a[i] 1 

end for 

- l;t < r;t • 

Si p[i]y[i] 


for {t 

e < 
u ■ 

p- 


w 


t + 1) do 


Si^H 
Si Wli]Y[i] 

Si Y[i] 

j exp(;3W[i]/ro) 

Si “[fiexp(/3rV[i]/W)^''F“ 


X uix + (1 — u!)x where x are tt, p'^, p^,u 
i* t— min (e*, 1.025e) 

Po t- p* + lO0,r(7r - TT*) + (fie log (e/i*) 

Po ■«- max(po,0) 
r •«— max{Q!r(^ — 7r),ro} 


V^£- 


0 


for (i 0; i < Af; z <— i + 1) do 
if a[i] == 1 then 

if £[i\ > —0W[i]Y[i] then 
f + max {£[i], 0} 



£ — min{f[z],0} 

m 

W[i]Y[i] 


$[ 1 ] ■<— min{^[z], 1/r} 

$[z] <— maxlihlz], —1/r} 

p+^pi(i-r$[z]) 
p- ^p“(i + r-i>[z]) 


0 Initial employment is 0.5 


l> binary variable: active (1) / inactive (0) firm 

t> Main loop 


0 Central Bank policy 


t> Firms npdate prices, productions and wages 


if y[z] < D[i\ then 
if 'P[z] > 0 then 

W[i] <— lF[z][l + 7 u,(l — rc|?[z])e random] 

W[i] ■<— min {W[z], (Flzj min [D[i], Y[i\\ + p'^ max {£[i], 0} + p^ min {£[i], 0})/y[z]} 

end if 

Y[i\ Tr- Y[i\ + min{77+(H[z] — y[z]),u*[z]} 
if p\i] < p then p[z] <r- p[z](l + 7p random) 

end if 

else if y[z] > D[i\ then 
if 'P[i\ < 0 then 

iy[z] •<— iy[z][l — 7 u,(l + r$[z])u random] 
end if 

y[z] max{0,y[z] — p_(D[z] — YJz])} 

if p\i] < p then p[z] t— p[z](l — 7p random) 
end if 
end if 

else if £[i] < —0iyjz]y[z] then 
a[z] 0 
V £[i] 

end if 
end if 
end for 
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Ei P[i]y[»l 

Ei yw 


/ = po + (1 - miE- 
d _ p’-e---v 
P S+£+ 

5^(1 + p‘')5 + Ei W\i]Y\i] 
c ■«- co[l + ac{n - p^)/7,r] 

Cb ^ cS 


Algorithm 2 MarkO (continued) 

l> Update u and p 


o Private bank sets interest rates 


[> Households decide the demand 


for (i 0; i < Ap; i <— i + 1) do 

nM J— Cbo-M exp(-Pp[i]/p) 

^ ‘ P[*] Ei “H exp(-Pp[ii/p) 

end for 

£+ c- 0 

for (i <— 0; i < Ap; z i + 1) do 
if a[i] == 1 then 

S <— S — p[i] min{y[i], -D[i]} 

P[z] p[i] minlyli], _D[z]} — iy[i]y[i] + max{£i[i], 0} + p^ min {5[i], 0} 
£[i] £[i] + P[z] 

if P[z] > 0 && f[z] > 0 then 

S ^ S' + 5£’[z] 

£\i\ £[{\ — 5£[i] 

end if 

max{f [i],0} 

end if 
end for 

7^ ^ 0 

for (i •<— 0; z < Ap; z z + 1) do 
if a[z] == 0 then 

if random < ip then 
y[z] •<— u random 

a[z] <r- 1 

P[z] <- p 
Ty[z] ■<— w 
£\i] ^ W[i\Y\i] 

TZ TZ + £{i] 

£'^ £^ -\- max {£[i], 0} 

end if 
end if 
end for 

for (z 0; z < Ap; z •<— z + 1) do 
if a[z] == 1 then 
if £[i\ >0 then 

£[i] ^ f [z] - n£[i]/£+ 

end if 
end if 
end for 
end for 


[> Inactive firms have no demand 

[> Accounting 


l> Pay dividends 


> Revivals 
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